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We present a conceptual design for an 
experiment to measure the neutron lifetime 
(-886 s) with an accuracy of 1 0^. The 
lifetime will be measured by observing 
the decay rate of a sample of ultracold 
neutrons (UCN) confined in vacuum in a 
magnetic trap. The UCN collaboration at 
Los Alamos National Laboratory has 
developed a prototype UCN source that 
is expected to produce a bottled UCN 
density of more than 100/cm^ [1]. The 
availability of such an intense source 
makes it possible to approach the 
measurement of the neutron lifetime in a 
new way. We argue below that it is 
possible to measure the neutron lifetime to 
10^ in a vacuum magnetic trap. The 
measurement involves no new technology 
beyond the expected UCN density. If 
even higher densities are available, the 
experiment can be made better and/or less 
expensive. We present the design and 
methodology for the measurement. The 
slow loss of neutrons that have stable 
orbits, but are not energetically trapped 



would produce a systematic uncertainty 
in the measurement. We discuss a new 
approach, chaotic cleanhig, to the 
elimination of quasi-neutrons from the 
trap by breaking the rotational symmetry 
of the quadrupole trap. The neutron orbits 
take on a chaotic character and mode 
mixing causes the neutrons on the 
quasi-bound orbits to leave the trap. 
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1. Introduction 

The beta-decay lifetime of the neutrons has a direct 
impact on cosmological models as a production of light 
elements during the Big Bang. The isotopic ratios 
measured in extragalactic gas clouds can be compared 
to Big Bang Nucleosynthesis calculations in order to 
place limits on the ratio of baryons to photons in the 
early universe. In these calculations of the expected 
"^He/^H ratio, the dominant uncertainty is the lifetime of 
the neutron [2,3]. 

Beta decay of the neutron is both the simplest 
nuclear beta decay and the simplest of the charged- 
current weak interactions in baryons. The weak inter- 
action parameters can be measured using neutron beta 



decay with fewer and simpler theoretical corrections 
than measurements using the beta decay of nuclei. The 
neutron beta decay rate is proportional to the quantity 
Sn + '^Sk where gy and gA are the semileptonic vector 
and axial- vector coupling constants. To extract the 
coupling constants from neutron beta decay measure- 
ments requires either an independent measurement of 
gv and gA or the ratio gjg^ = X, 

The other important goal of measuring the neutron 
lifetime with improved accuracy is to test unitarity of 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix. For 
the unitarity the sum of the squares of the first row 
CKM matrix elements must be one 



+ KJ'+K. 



1. 
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Experimental constraints on the values of X and | V^^ 
are coming from nuclear 0^ - 0^ beta decays [4], CKM 
unitarity (assumes conserved vector current, CKM 
unitarity and values of \V^ and \V^ from [3]), meas- 
urements of the neutron beta decay correlation A [5] 

A = -1.2735 ±0.0021 
|FJ = 0.9756 ±0.0005. 

In order to improve the determination of | V^^ beyond 
the precision of obtained from the above systems, one 
can study the neutron beta decay. An advantage of the 
study of the neutron system is the relatively stronger 
understanding of the necessary corrections. The nucleus 
dependent radiative correction 5r has been calculated to 
the 10"^ level for the neutron [6], and there is no 
Coulumb correction. Thus, the limiting theoretical un- 
certainty on \V^ determined from neutron decay 
is simply that due to the inner radiative correction, 
zi\(± 0.0004) [4]. The overall radiative corrections to 
the neutron decay rate have been calculated at the 1 0"* 
level [6,7]. At this level, the neutron beta decay rate can 
be written as [7] 

r= 0.1897|FJ'(l+3A')(l+0.0739± 0.0008) X 10"' s-\ 

The CKM matrix element | V^^ can thus be determined 
from the neutron system by measuring the neutron life- 
time (1//") and X, The ratio X is most precisely deter- 
mined from measurements of the electron-neutron 
spin-asymmetry coefficient A, Currently the limiting 



uncertainty on \V^^ measured with neutrons comes 
from A. As measurements of the electron-spin correla- 
tion, A, and the electron-neutrino correlation, a, are 
improved [8], towards the goals of proposals to meas- 
ure them [9], the neutron lifetime will limit the accura- 
cy of CKM unitarity tests. A measurement of the neutron 
lifetime with an accuracy of 10"* will be more than 
sufficient for astrophysics and for CKM unitarity tests. 
The accuracy of the measurement will challenge the 
theory of inner radiative corrections. 

2. Vacuum Quadrupole Trap 

Figure 1 shows schematically the proposed quadru- 
pole trap geometry and Fig. 2 shows the trapping field. 

The UCNs have velocities of up to 5 m/s and the 
strength of the trapping field is 2.2 T. The projection of 
the neutron magnetic moment on the magnetic field 
direction is an adiabatic invariant. For fields of the 
order of tesla, depolarization is negligible, 10"^^ in 10"^ s. 
The above potential has minimum of |5| and one spin 
state is trapped and the other expelled. Lowering 
current in loop 1 fills the trap. The material guide deliv- 
ers UCNs into the trap. In a few seconds the neutron 
population in the trap comes into equilibrium with the 
flux from the guide and the current is increased to close 
the trap. The neutrons move in the vacuum of the warm 
bore of the superconducting quadrupole magnet. Once 
the trap is closed there are no losses from the trap other 
than neutron (3 decay and the possible loss of quasi- 
trapped neutrons discussed below. 




Guide 
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Black Absorber 
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Fig. 1. Proposed trap geometry. 



362 



Volume 110, Number 4, July- August 2005 

Journal of Research of the National Institute of Standards and Technology 



trapping field vs , r 



trapping field vs » z 
.4 




-0^1 -Q.OS 0,05 0,1 -0.75 -0.5 -D.25 0.25 0.5 0.75 

Fig. 2. Trapping field, \B \ in tesla, versus the axial coordinate, z in meters, and the radial coordinate, r in meters. 
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Fig. 3. The surface of |5 1 = B^^ or U= U^^ for the quadrupole trap and the contours of \B \ at the center of the trap. The neutrons escape 
the ridge at the top of the trap — the saddle point on the contour plot — or at the cusps at the ends of the trap. 



along 



The above trap has been designed to have a low 
symmetry for two reasons. First, a low symmetry 
reduces the probability of quasi-trapped orbits by 
inducing mode mixing of the neutrons so that neutrons 
reach the boundary defined by their kinetic and poten- 
tial energy Second, the approximately square trap cross 
section helps to match scintillation detectors to the trap 
shape. Neutron decays are observed by detecting the up 
to 0.78 MeV p particles from neutron decay. 



The trap geometry is chosen to facilitate the detec- 
tion p particles. When a neutron decays in the trap, the 
emitted p spirals around field lines. Since the field lines 
end on the quadrupole pole faces the ps are guided to 
the poles. Cosmic-ray events are vetoed by the veto 
scintillators shown in Fig. 4. 



363 



Volume 110, Number 4, July- August 2005 

Journal of Research of the National Institute of Standards and Technology 



Veto y'^osmic 




Field Line 
p Detector 



|i3| from all 



0,4 



0,2 




Fig. 4. Detector scheme. Four comer-shaped scintillation (blue) detectors cover the four comers of the trap. Veto scintillators (green) to reduce 
cosmic ray backgrounds surround the detectors. We estimate that the veto inefficiency will be less than 10"^ and the cosmic/signal ratio will be 
1/10. The ps spiral around field lines and are guided to the detectors. About 20 % of the ps that strike a detector will be reflected. These back- 
scattered ps again follow field lines and strike the other side of the same detector. The four veto detectors and the other three scintillators form an 
effective veto of cosmic rays and other backgrounds. On the right are shown the contours of constant yield from the trap. The density of trapped 
neutrons is high in the center of the trap, but the magnetic pinch effect prevents ps from these neutrons from reaching the detector. There are few 
neutrons trapped near the outer edges of the trap. Most detected ps come from the mid-field portion of the trap. About 30 % of the trapped ps are 
detected. The horizontal and vertical units on the right-hand contour plot are one graph unit = 0.1 m. 



3, Elimination of Quasi-Bound Neutrons 

We have studied two approaches to the elimination 
of quasi-bound orbits in a two-dimensional (x~y) trap. 
Both approaches work by breaking the symmetry of the 
trap in a way that causes the neutron orbits to take on 
a chaotic character. Quasi trapping occurs when the 
energy of the neutron is shared between modes. The 
kinetic energy of the neutron never becomes zero and 
the neutron never reaches the surface U= U^^^. In the 
first approach, shown in Fig. 5, we add a line current 
in the z direction to the two-dimensional quadrupole 
trapping potential. The axial symmetry of the trap is 
broken, angular momentum is not conserved, and the 
orbits take on a chaotic character. 

In the second approach, we insert a ^^Ni mirror that 
reflects neutrons. The mirror may be inserted like a 
knife so that the energy added to the neutrons that 
collide with the mirror while it is being inserted is 
small. Between collisions, neutrons move in orbits 
having constant angular momentum. When a neutron 
collides with the mirror, its angular momentum changes. 



If the angular momentum falls below some critical 
value, the neutron crosses the circle of r+ V=U^^^ and 
is removed from the trap. 

Quasi-bound orbits may be eliminated by adiabati- 
cally reducing the strength of the trapping field, how- 
ever during this procedure a large fraction of energeti- 
cally trapped neutrons are lost along with the quasi- 
trapped neutrons. In the above example the field 
strength must be reduced by a factor of 8/27, and 71 % 
of the energetically trapped neutrons would be lost. On 
the other hand, the time required for chaotic cleaning of 
the quasi-bound orbits increases as the ratio of the 
neutron energy to U,^ approaches unity. The best 
approach may be to first clean the trap chaotically for a 
few seconds and eliminate orbits with t/> (1 + s)U^^ 
where £^0.01. Then the quasi-trapped orbits can be 
cleaned by a small field reduction (B -> 0.97 B^^). 
Furthermore, it may require a time larger than the 
neutron lifetime to substantially lower the trap field 
strength. The smaller field strength reduction needed 
for chaotic cleaning is an important advantage of 
chaotic cleaning. 
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Fig. 5. The upper row shows a symmetric quadrupole potential. The left plot shows trapping field, \B\ in tesla versus x andy. The spatial dimen- 
sions on all plots are one graph unit = 0. 1 m. A quasi-trapped orbit is shown in the right plot. The lower row is for the same quadrupole with a line 
current added to break the symmetry. The depth of the traps, the energy of the neutron, and the initial position and velocities are the same. The 
orbit crosses curve of T+V= t/j^^ (between light blue and light green on the contour plot and the red line on the orbit plots) for the lower 
chaotic trap but not for the upper symmetric trap. 



A problem that must be addressed is the activation of 
the guide and other objects when the trap is filled. If the 
act of filling the trap produces radioactivity with a life- 
time comparable to the neutron lifetime, and the decays 
are detected, a systematic uncertainty in the measured 
lifetime will result. The trap and guide system will be 
designed to minimize activation, but the background 
must be measured in situ. The following procedure will 
be used to mitigate this effect. First fill the trap as if to 
begin a measurement cycle, but then lower loop 2 and 



allow the trapped neutrons to escape into the black 
absorber. Then restore loop 2 and measure the back- 
ground. 

We estimate that 8 x 1 0"^ p decays are detected per 
fill in a 27 L trap. The neutron lifetime could be meas- 
ured with a statistical uncertainty of 10^ in 34 days 
with an UCN density of 100/cml If a higher UCN 
density were available, the measuring time could be 
shorter or the high density could be used to reduce the 
size of the trap, the field strength, or in other ways. 
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fit and data vs bin numbeir 
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Fig. 6. Upper two plots show trajectories in a two-dimensional quadrupole trap. The units for the horizontal and vertical axes are meters. The left 
figure shows the motion of a neutron with energy 1.02 times t/niax- The neutron never reaches the (blue) circle T+ V= t/niax- When a (blue) 
mirror is inserted as shown in the upper right figure, the angular momentum is modified by coUisions, and the neutron orbit crosses the circle after 
a few periods and is removed from the trap. The lower graph shows the distribution of removal times. The size of the time bins for the x-axis is 
0.3 s. The approximately exponential part of the distribution has a time constant of 0.54 s. The time required for the quasi-trapped population of 
neutrons to be reduced to 10"^ of its initial value is 6.2 s. 
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